Understanding self-assembled nanosphere patterns 
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Patterns generated by a colloidal suspension of nanospheres drying on a frictional substrate are 
studied by experiments and computer simulations. The obtained two-dimensional self-assembled 
structures are commonly used for nanosphere lithography. A spring-block stick-slip model is in- 
troduced for simulating the phenomenon and the influence of several controllable parameters on 
the final structure is investigated. The model successfully reproduces the experimentally observed 
patterns and the dynamics leading to pattern formation is revealed. 



The so-called bottom-up approach for the fabrication 
of nanostructures starting from stable building blocks 
such as molecules or nanoparticles has become an in- 
creasingly popular topic in nanoscience and nanotechnol- 
ogy. Thanks to the efforts of nanochemists, during the 
past decades various nanoparticles of polystyrene, silica, 
noble-metal and semiconductor, nearly monodisperse in 
terms of size, shape, internal structure, and surface chem- 
istry, can be produced through a reliable, standard manu- 
facturing process. Using these nanoparticles as building 
blocks, the synthesis of long-range-ordered monolayers 
and films of colloidal nanocrystals has been in partic- 
ular focus. The revolutionary development of photonic 
crystals triggered efforts to get innovative methods for 
crystallizingpolystyrene colloids and creating new crystal 
structures yj, |2| • The use of two-dimensional (2D) self- 
assembled array of nanometer-sized polystyrene spheres 
as deposition mask is known as NanoSphere Lithogra- 
phy (NSL) Q . The homogeneous arrays of nanoparticles 
produced using NSL are potentially useful in studies of 
size-dependent optical, magnetic, catalytic and electri- 
cal transport properties of materials 0, S 0, E|- NSL is 
now recognized as a powerful fabrication technique to in- 
expensively produce nanoparticle arrays with controlled 
shape, size and interparticle spacing. A fundamental 
goal for further progress in NSL is the development of 
experimental protocols to control the interactions, and 
thereby the ordering of nanoparticles on solid substrates 
0, @- However, it is more and more clear that due 
to the rich physics and chemistry underlying the forma- 
tion of nanoparticle arrays from colloidal suspensions, the 
likelihood of structures other than close-packed networks 
forming during solvent evaporation is very high 0, llf]j | . 
Therefore a major motivation for theoretical research in 
this field remains the challenge to understand how or- 
dered or complex structures form spontaneously by self- 
assembly, and how such processes can be controlled in 
order to prepare structures with a pre-determined geom- 
etry 1 1 1 1 ] . The present study intends to contribute in this 
sense by proposing a model that can be easily studied 
through computer simulations and it is able to qualita- 
tively reproduce the wide variety of observed patterns. 
We focus on an experimentally simple case, when 2D self- 
assembled arrays of nanometer-sized polystyrene spheres 






FIG. 1: Characteristic two-dimensional patterns obtained 
from a colloidal suspension of 220 nm diameter nanospheres 
drying on a silica substrate. 



will form from a colloidal suspension which is drying on 
a substrate. Some characteristic patterns obtained as a 
result of this phenomenon are present in Figure ^ 

The model is rather similar with the spring-block stick- 
slip model successfully used for describing fragmentation 
structures obtained in drying granular materials in con- 
tact with a frictional substrate 0, 0] . The new feature 
of the present model is that a predefined lattice is not 
considered anymore. Proceeding in such manner the fi- 
nal lattice structure appears as a result of the underlying 
interactions governing the dynamics of the system and 
the existing geometrical constraints. 

The model is two-dimensional; its main elements are 
blocks which can move on a frictional substrate and 
springs connecting these blocks (Figure |2Jl). Disks, 
all with the same radius R, model nanospheres, while 
springs which are characterized by their variable length 
and similar spring constants k, represents water bridges 
between them. The length of a spring is defined as the 
distance between perimeters of connected disks. There 
is also a Lenard- Jones type interaction- force Fj , between 
each pair of disk. This is characterized by a strong, al- 
most hard-core type repulsion which forbids disks to in- 
terpenetrate each other and by a weak attractive type 
force, accounting for the electric Van der Waals type in- 



2 



2R 



2R 

K 



If. 



Lena rd- Jones type p otential 



F. 



Id) 




FIG. 2: Sketch of the spring-block stick-slip model. Fig- 
ure (a) and (b) shows the distance dependence of the ten- 
sion Fk in the spring and the Van-der Waals type force Fj 
between nanospheres, respectively. Figure (c) illustrates the 
friction force Ffric between nanospheres and substrate, as a 
function of total force Ft acting on nanospheres. Figure (d) 
presents the main elements of the model, and Figure (e) is an 
example for a simplified spring-block monolayer system. 



teraction between nanospheres (Figure |2p)- The friction 
between disks and surface equilibrates a net force less 
than Ff (Figure Et) ■ Whenever the total force acting on 
a disk exceeds Ff, the disk slips with an over-damped 
motion. The tension in each spring is proportional with 
the length of the spring (Fk = k ■ I), and has a breaking 
threshold F b (Figure UK)- 

Initially disks are randomly distributed and connected 
by a network of springs (Figure We put springs 

between those spheres, for which their center can be con- 
nected without intersecting another sphere (this condi- 
tion will be referred later as the geometric condition). 
An initially pre-stressed spring-block network is thus con- 
structed. During each simulation step the spring con- 
stant is fixed and the system relaxes to an equilibrium 
configuration where the tension in each existing spring is 
lower than the breaking threshold Fb, and the total net 
force acting on each disk is lower in magnitude than the 
slipping threshold Ff . This relaxation is realized through 
several relaxation steps. The time length dt for each re- 
laxation step is taken as unity (dt — 1). (Since we are 
interested mostly in the final, stable structure, connec- 
tion with real time is not crucial for us, and we can define 
time units in an arbitrary manner.) Considering a clas- 
sical molecular dynamics simulation for relaxation would 
be very time-consuming. Following the method used for 
simulating drying processes in granular media |l2l llflj , 
we choose thus a simplified dynamics, where the connec- 
tion with real time is lost, but the relaxation remains 
realistic: 



(1) For all springs the tension 



is compared with 
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from the system. 

(2) The total forces F\ 



> Fb, the spring is broken and taken away 



Ff ) acting on 



disks are calculated (the sum is over all the other disks p, 
dip is 1 if the disks are connected by a spring and oth- 
erwise, the subscripts k and j denotes elastic forces from 
springs and Van der Waals type forces between disks, 
respectively) . 

(3) Each disk is analyzed. If the magnitude of the 
total force F[ acting on a disk is bigger than Ff, then 
the disk will slip with an over-damped motion governed 

-* i 

by viscosity r\, and its position will be changed by: dr = 
Fldt/rj. The repulsive part of the Lenard- Jones potential 
forbids the spheres to slide on each other and the presence 
of viscosity eliminates unrealistic oscillations. 

(4) During the motion of a disk it can happen that 
another spring is intersected. This intersected spring will 
brake and it will be taken away from the system. 

(5) After all disks have been visited and their possi- 
ble motions done, the springs that fulfill the considered 
geometrical condition and for which the tension is lower 
than the breaking threshold are redone. By this the re- 
arrangement of water between nanospheres is modelled. 

This concludes one relaxation step. The relaxation is 
continued until a relaxation step is finished without hav- 
ing any spring breaking or disk slipping event. Since in 
our algorithm it takes a very long time to achieve a per- 
fect relaxation, we introduce a tolerance, and assume the 
relaxation completed when the maximal slip is smaller 
than this tolerance. After the relaxation is done, we pro- 
ceed to the next simulation step and increase all spring 
constants by an amount dk. This models that water 
bridges are getting thinner and the capillary forces in- 
crease the tension in them as the water evaporates. The 
system is relaxed for the new spring-constant value, and 
the spring constant is increased again, until all springs 
are broken or a stable limiting configuration is reached. 

The above dynamics can be easily implemented on 
computer and relatively big systems with over 10000 of 
disks can be simulated in reasonable computational time. 
Several types of boundary conditions can be considered. 
One possibility is the use of free boundary condition 
which can be realized in a simple manner by positioning 
initially the disks inside a circle to minimize the effect of 
edges. Another possibility is to consider fixed boundary 
conditions. This can be realized for example by posi- 
tioning again the disks inside a circle and considering a 
chain of fixed disks on the chosen circle. These fixed disks 
are then connected with geometrically allowed springs to 
other disks. One can also consider periodic boundary 
condition and position initially the disks inside a rectan- 
gle. It is quite easy to check that the considered bound- 
ary condition will influence the final stable structure only 
in the vicinity of boundaries. In the bulk, far from the 
boundaries, the obtained structures are rather indepen- 
dent of the imposed boundary conditions. Simulation re- 
sults presented here are taken from a small central part 
of systems with N > 10000 disks using fixed boundary 
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conditions. 

The model, as described above, has several parameters: 
the value of the static friction Ft, the spring breaking 
threshold value Fb, the initial value of spring constants 
ki n i, the spring constant increasing step dk, the viscos- 
ity 77, the parameters of the Lenard-Jones potential, the 
radius of disks R, and the initial density of nanosphcrcs 
q = S/ (NttR 2 ) (where S is the simulation area). By 
direct simulations one can verify that if small enough dk 
is chosen, in the quasi-static limit the final structures are 
not sensitively influenced by this parameter. The model 
will only work for viscosity values chosen between rea- 
sonable limits, and for these viscosity values the final 
patterns are rather similar. (Choosing a too small vis- 
cosity will result in unrealistic oscillations of disks, while 
a too high value will make the disk slip too small and 
increase considerably the simulation time). It is desir- 
able to choose the value of fci„i small, to simulate the 
whole stress building process in the system. The Lenard- 
Jones potential is fixed, so that the repulsive part gives 
a strong hard-core repulsion for distances smaller than 
2R, and a small attractive contribution to the net force 
between pairs of disks for distances bigger than 2R. The 
radius of disks were considered as unity (i? = 1), defin- 
ing the unit length in the system. We remain thus with 
three main parameters governing the generated patterns: 
i 7 },, Ff and g. The influence of these three parameters 
on the final structure was investigated and the obtained 
structures were visually compared with experimentally 
obtained ones. It is worth mentioning here, that for a 
successful simulation the values of g and Fb, has to be 
correlated. For small g values the value of Fb has to be 
chosen high enough, otherwise the majority of springs 
would break in the very first simulation step, and no 
clusterization would be observable. 

The experimental samples were prepared using the 
drop-coat method 0- 0|. Of critical importance for 
nanosphere ordering is the chemical treatment of the 
glass substrate to render the surface hydrophilic and im- 
prove its wettability. This was achieved by etching the 
substrate in a solution of sulphuric acid/hydrogen perox- 
ide (3:1) for a period of 3 hours. The substrates were then 
washed in a copious amount of deionized water, immersed 
in a solution of ultra pure water/ ammonium hydrox- 
ide/hydrogen peroxide (5:1:1) for 2 hours and sonicated 
in an ultrasonic bath for 1 hour. Finally, the substrates 
were thoroughly rinsed in ultra pure water and stored 
under ultra pure water. Polystyrene nanospheres of 220 
nm diameter, exhibiting negatively charged carboxyl- 
terminated surface with a strongly hydrophobic nature, 
were supplied as monodispersed suspensions in deionised 
water (wt 4%). The original suspension of polystyrene 
nanospheres was diluted by 10 and a volume of 100 1 di- 
luted solution was evenly spread on the pre-treated sub- 
strates. Finally, samples were dried in an oven at 65 °C 
for 45 minutes. As the water evaporates the nanospheres 
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FIG. 3: Characteristic simulation results for the influence of 
nanosphere density on final patterns. 

self- assemble into close-packed monolayer arrays exhibit- 
ing many fracture lines, dislocations and other type of 
defects (Figure^]). The microstructure of the sample was 
studied by scanning electron microscopy (SEM) using a 
JEOL JSM 5600 LV electronic microscope. By surveying 
the substrate, surface regions with different morpholo- 
gies were found even within the same sample. These 
qualitatively different structures will arise mainly due to 
the non-homogenous nanosphere concentration in differ- 
ent parts of the sample. 

On Figure |21 we show computer simulation results for 
the influence of nanosphere density on the final struc- 
ture. For high densities, in agreement with experimental 
results, we obtain close-packed domains (triangular lat- 
tice structure), separated by realistic fracture lines and 
defects. As the nanosphere density decreases we repro- 
duce the experimentally observed isolated and elongated 
islands. 

Similarly with the results obtained by Leung et al. 01 
we observed that the topology of the final structure is 
strongly affected by the values of Ff and Ft,. Increasing 
the value of Ff will result in denser and thinner fracture 
line structure with smaller close-packed domains (top line 
of Figure 0J|. Decreasing the value of Fj, results in a 
similar effect (bottom line of Figure 0J. Experimentally 
one could control Ff by changing the substrate and Fj, 
by changing the evaporating liquid. 

The dynamics leading to formation of final structures 
can be also easily revealed through simulations (to do 
that experimentally is quite a difficult task). Some 
movies in this sense are p resented on the home-page ded- 
icated to this study Jljj. The formation of close packed 



4 



i' 



- i < •••57 




F b = 0.02 



F b = 0.05 



F b = 0.1 



FIG. 4: Simulation results investigating the influence of pa- 
rameters Ff and Ft on the final patterns. On the top line 
Ff is changing (Fb = 0.05), while on the bottom line the var- 
ied parameter is Fb (Ff — 0.005). For all these simulations 
the other parameters were chosen as: kmi — 0.1, dk — 0.002, 
g = 0.749, 77 = 250. 




FIG. 5: Experimental (top) and simulation results (bottom) 
for the square lattice crystallization phase. 



(triangular lattice) domains separated by characteristic 
fracture lines is realized through a series of stick-slip 
avalanches, as the value of the spring constant is in- 
creased. The initial nanosphere configuration is quickly 
lost, and due to long-range elastic forces an interme- 
diate square lattice multi-domain structure is formed 
(Figure bottom line). This intermediate (metastable) 
structure is lost as k increases further (long springs are 
broken) and slowly the triangular symmetry is selected. 



During the relaxation of stress and non-homogeneity in 
the system, fracture lines are nucleated and propagated. 
This process will lead to the final characteristic pattern. 

The intermediate square lattice structure can be also 
revealed by experiments, either by analyzing samples 
that are not completely dried, or by considering a quick, 
non quasi-static drying procedure. It is believed |l5j that 
in case of non quasi-static drying the nanosphere system 
is blocked in this intermediate, metastable phase. Exper- 
imentally, we have also obtained samples that show this 
intermediate square-lattice crystallization phase (Fig- 
ure [SJ top line). The success in reproducing this square- 
lattice crystallization phase and in general the qualita- 
tive agreement between simulation results (Figure I2l.'{ll 
and the experimentally observed structures (Figure ^) 
encourages us to believe that the introduced model is a 
good one and the picture we got for the underlying dy- 
namics is correct. 

We have introduced thus a successful model for de- 
scribing monolayer patterns that are formed on a surface 
after quasi-static drying of a colloidal nanosphere suspen- 
sion. By using this model one can get a first picture on 
the underlying nano-scale dynamics and can analyze the 
influence of the experimentally adjustable parameters on 
these practically important structures. 
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